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TECHNICAL MEMORANDUM 


CLOSED-FORM SOLUTIONS FOR ATMOSPHERIC FLIGHT 
WITH APPLICATIONS TO SHUTTLE GUIDANCE 

INTRODUCTION 

The numerical integration of any set of differential equations for the 
atmospheric motion of a rocket-powered vehicle is usually a time-consuming 
procedure even when a high-speed digital computer is used. For this reason, 
approximate closed-form solutions, which take much less time to evaluate, are 
often used to represent the atmospheric motion of a rocket-powered vehicle. A 
set of closed-form solutions for reentry and atmospheric ascent is developed 
in this report. For both cases accurate closed-form solutions for the three- 
dimensional Cartesian coordinates of the vehicle’s position and velocity vector 
are obtained. Also, the partial derivative transition matrices of the final 
state of the vehicle with respect to the initial state are obtained as closed-form 
expressions. In addition, the closed-form solutions for the state and partial 
derivative transition matrices are used to develop a variation-of-parameters 
technique which can provide rapid numerical integration of the original 
equations of motion. 


DEVELOPMENT OF CLOSED-FORM SOLUTIONS 

The starting point for all the work in this report will be the equations of 
motion for atmospheric flight developed in Reference 1. The techniques devel- 
oped could be applied to any other equations of motion, but a particular example 
is useful in illustrating the ideas. Thus, the equations of motion derived in 
Reference i are chosen as an example for this report. The derivations of these 
equations will not be repeated here, therefore, the reader may want to 
examine iReference 1 before starting the next section. 

The equations of motion for the center of gravity of a space vehicle in 
an inertial three-dimensional Cartesian coordinate system as given in Refer- 
ence 1 are first-order ordinary nonlinear differential equations for the position 
vector (x of dimension 3) and the velocity vector (v of dimension 3); i. e., 

x = v 


F_ 

m 



T, + D 


GMT 

R 3 


v 


m 
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where 


F = F + A (P 0 - P) , 

D C 

m = mo - m(t - t 0 ) , 

X “I P A r °L a [ l 7 r i! (i|r) - ( T r ' m)\ ] • 

B = -(? p A r) ['V ( C A + 2ilC lJ ' ' ifl)] 7 r ’ 

and 

V = 7- 5) Xx-W 
r 

It should be noted that p is the three dimensional control vector which must be 
determined or specified as a function of time, and the quantities F , A , P 0 , 

__ o 6 

m 0 , m, A^, oj, and GM are constants defined in Reference i. The quantities 

p, P, c , c , and rj are (in general) specified functions of x and v. The 
A L 

G? 

exact form for these functions is also given in Reference 1. With the pre- 
ceding statements, it can be seen that the right-hand side of the v equation is 
a function of t, x, and Thus, a given time function for the control vector 
p and initial conditions x 0 and v 0 at some time t 0 allows numerical integration 
to be used to obtain values of x and v at any t > t 0 . Usually the time function 
for control vector p is determined by optimization theory to maximize or 
minimize a quantity J defined by 

J = <p (x f , v f , y + cq f 

and to satisfy a vector function (of dimension < 6) called the physical bound- 
ary condition and written in the following general form : 

F(x f , v f , t f ) ■ 0 

In the above expression for J, c is an arbitrary positive or negative weighting 
factor for q^; and q^ is a new variable defined by another first-order nonlinear 

ordinary differential equation with the following form: 
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q = f(x, v, p, t) . 


Then, for a given control vector time function and initial conditions x 0 , v 0 , 
q 0 , at some time t^, numerical integration can be used to obtain x^, v , and 

at the final time t^. Thus, it can be seen that the value of J depends on the 

control vector time function and the initial conditions x 0 , v 0 , and q 0 at t 0 . 

Since the initial conditions are usually assumed to be specified by the nature of 
the problem, the only thing left that affects the value of J is the form of the 
control vector time function p(t), where to ^ t ^ t • 

In order to use optimization theory to determine the control vector 
time function p(t), the time functions A (a three-dimensional vector), u 
(another three-dimensional vector), and y (a scalar) are introduced. These 
time functions are usually called Lagrangian multipliers, or simply multipliers, 
and are defined by their differential equations. To obtain the differential 
equations for A, u, and y, define 

-T • _T • 

H=Ax+. uv+yq, 

where the superscript T denotes the ordinary matrix transpose. 

Then 

iT_> 8H 
X ~ " 3x ’ 

j_T 3 H 
u = - -rr . 



The above expressions allow the differential equations for A, u, and y to be 
obtained by performing the indicated partial differentiations of H with x and v 
replaced by their right-hand sides in the expression for H. Then values for 


”X., n, and y at any t > t 0 can be obtained by numerical integration of the 

-s-T -lT 

differential equations A , u , and y along with a simultaneous numerical 
integration of the differential equations X, v, and q when a control vector 
time function p(t) , where t 0 ^ t — t^, is specified and the initial conditions 

* 0 , 3o. Yo> *0* 7 0 » flo at to are specified. 
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Now the results of optimization theory can be used to state that any 
control vector time function which minimizes J must maximize at every t, 
where t 0 — t — t^, the quantity H, and any control vector time function which 

maximizes J must minimize at every t, where t 0 < t ^ t^, the quantity H. 

In addition, any control vector time function which minimizes or maximizes J 
and satisfies the physical boundary conditions [denoted by F(x, v, t ) = 0] 

must maximize or minimize H as stated previously and must be produced by 
initial conditions Uq, and y 0 that satisfy the following conditions: 



H = 0 . 

f 


In the above expressions, ’fr is another vector of constant multipliers that 
must be introduced and has the same dimension as F(l?, v, t f ). If a control 

vector time function and the multipliers A. 0 , Uq, y 0 , and 'k can be found that 
satisfy the preceding necessary conditions, the optimization boundary valve 
problem has been solved. Reference 1 discusses this problem in more detail 
and indicates a few techniques for solving the problem when numerical integra- 
tion is used on all the ordinary nonlinear differential equations involved. The 
purpose of this report is to make simplifications in the ordinary nonlinear 
differential equations involved so that closed-form solutions can be used 
(instead of numerical integration) to solve the optimization boundary value 
problem. This work will be done separately for the ascent case and the reentry 
case. 


Ascent Case 

For the ascent case, let q = m and J = -q^ so that it is desired to 

determine the control vector time function 'p(t) to minimize J, thus maximiz- 
ing mass. Then the equations of motion become 

x = v 


4 



• F / p \ L + D GMx 

7= qtfW ~ *r • 

q ~ ~m 

When the expressions for X. and D are substituted into the T equation, the 
following expression is obtained: 



This expression can be rewritten as follows: 
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Then 


g = - 


GMx 

R S " 


v = 
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so that the expression for H can be written easily as 

H»* T V+n T {(±) [A] (i)S+gj -ym • 

Then 



Since J is to be minimized, p must be chosen to maximize H at every 
t, where t ^ t ^ t^. By examining the expression for H, it can be seen that 

_ T — — 

p = [A] u is sufficient to assure a maximization of H with respect to p at 

every t such that t^ — t — t^. This result for p can be substituted back into 

all the differential equations, and the formulation of the optimization boundary 
value problem is practically complete. The entire system of differential 
equations that must be integrated simultaneously can be written as: 



q = -m , 


6 




/ 1 \ / / U T [AI \ 9([AI X u) , -T 3b 
W IW[A][A] T u/. 8 * " 8 * 


-Tag 
u ttS- 
a x 


* 



—T 

u [A] 


Vi 1, 


(AHA) 


-) 


a([ a] T u ) 
a v 



» 



(A] 


/ U] T u \ , 
' n/ u T [ A] [A] u' 


Given initial conditions 3^, V 0 , q 0 , X 0 , TT 0 , and y 0 at t 0 , this system 
can be integrated numerically; but a closed-form solution is to be derived here 
that will approximate with some degree of effectiveness the solution obtained by 
numerical integration. The first step in obtaining a closed-form solution for 
the entire system is to write X{t) and TT(t) as a Taylor Series expansion about 
the initial conditions; i. e., 

T(t) = T 0 + "^(t - t 0 ) + — Afl(t - t 0 ) 2 + ... 

u(t) = Uq + u 0 (t - to) + i %(t - to) 2 + . . . . 

±. ^ ii. 

The differential equations for X and H can be used to evaluate u^, X*, 

Tifl, etc. Also, the quantities [A] , b, and g are expanded in a Taylor Series 
about the initial conditions; i. e. , 

[A] = [ A 0 ] + [A 0 ](t - t 0 ) + j [A 0 ](t - t,,) 2 + ... , 

b = b 0 + t 0 (t - tg) + — ^j(t - tg) 2 + ... , 

g = go + £<>(* - to) + J fo(t - t 0 ) 2 + ... . 

The quantities [A 0 ], b 0 , g 0 , [A 0 ], T5 0 , g 0 , etc., are determined by differen- 
tiating the expressions for [A] , 15, and g with respect to t and evaluating 
them at i^. Then the Taylor Series expansions for [A] , b, gX, and u can 
be truncated at some power of (t - to) and substituted into the expressions for 
T and y to give expressions for T and y whose right-hand sides are functions 
of only the independent variable time. Then numerical quadrature techniques 
can be used to find the solutions v(t), y(t), and lf(t), where t 0 ^ t — t^, 
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much more rapidly than could be done by simultaneous numerical integration of 
the original system. The accuracy of this approach is limited only by the 
accuracy of the truncated Taylor Series expansions over the time interval of 
interest. To go further and obtain actual closed-form solutions requires a 
few more approximations. First X(t) and u(t) are assumed to be represented 
adequately by using only the linear terms in the Taylor Series expansions. 

Then [A] must be assumed to be constant. With these assumptions, the 
expressions for v and y become 


r "i Tr -i- 
-A n j Lu<i + up(t - t n )J 


, jv . * iTr^ , 

^ J u? + (t - to) [A 0 ] [A 0 ] T [uo r + ^'(t - to)] 


+ ^-[bp + b 0 (t - to) + ^ b 0 (t - t 0 ) 2 + ...] 


+ [ go + go(t - to) + \ go(t - t§) + . . . ] , 


1 

y - — 5 

q 



u 0 T + u 0 T (t - t 0 )J [a 0 ] [a 0 ] T [u 0 T + u„ T (t - t 0 )J 
[u? + u?(t - t 0 )J [Fo + t 0 (t - t 0 ) + | F 0 (t - t 0 ) 2 + . . .] | 


Since q = -m, then q = q 0 - m (t - tj so that = — r- and (t - t 0 ) 

dq m 

= -X (q - q 0 ). With these substitutions, and 4^ can be written 
m dq dq 

follows: 


as 


S--il(i) 


[ A o] [A.] [% - - <to)] 


J - §-(q - % ) [ A o] M T [^I - ^(q - % )] 


+ I 


b 0 - ^r- (q - q 0 ) + \ ^Mq - q 0 ) 2 - ... 

m 2 m 


ql_ m 


+ go - -§-(q - qo) + |§r(q - qo ) 2 - ••• j * 


m 
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+ -d&M - «o> 2 - 


Now these equations can be integrated fairly easily in closed form. To illus- 
trate this closed-form solution, only the constant terms in the Taylor Series 
expansion for 15 and "g will be used, although the integration for the inclusion 
of any number of terms in the series for b and g is easily added. First, note 
that q = m and q 0 = m 0 so that the following quantities can be defined: 


- [Ao ] l .'*oJ u o + Uq - f Aq 1 [ A# 3 TTq - • m o ^ 2 , 


= -[Ao 1 

[Ao3 T u 0 | 

» 

r_T , 

= [uo + 

*r®] 

[ Ao ] [Ao ] T Ju 0 + u 0 (-^) ] 

r_ T 

= [uo 

* 

] d t = u 0 T d t - de , 

f-T 
= -:[uo 

^ iT /mj\ 
+ U ° \m ) 

j [ Ao 1 1 Aq ] 

-m 

•T 

uq [ Aq ] Ac 

) T [uo + %(^)] 


= - 3 < r “»(s) - (i) d ’< = -(I ' ”» T d * 

dg = 1A 0 1[A 0 ] u 0 =-( s )u, dj , 
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d io = [uo T + £ (^)] b o . 

d n = -(^)^ T b ° * 


become 


With these definitions, the expressions for ^ and ^ 

dq dm dq dm 


dv _ / 1 \ I /_1\ dj + d^m 

to • -\Z) \ Km) j ar + ^ m + d>rf 


1\- . - 
go 


(3* 


dy 

dm 


= ‘(^') (^) { ^ d 7 + d8 m + dgm 2 + die + dum J 


dv dv 

The above expressions for ~~ and -7— can now be used along with a table of 

dm dm 

integrals to obtain expressions for V(m) and y(m) at any m > m D . If v(t) 
and y(t) are desired, then the relation m = m 0 - m (t - t 0 ) can be used to 
calculate the m corresponding to any t > t 0 . To see how the expressions for 
v(m) and y(m) are obtained with the aid of a table of integrals, define 


X = d 7 + dgm + dgm 2 


Then 


-Id, 

f m 1 

C dm 

+ (L 

■ m 

r 

dm 

+ b 

0 

’ m 

f 

dm 

m | 1 

d nr/lT 
Lmo J 


j 

-™o 

7T 

J 

L^o 

m 


r m 


go 


/ dm 
m 0 


1 f r m 

y(m) = y 0 - ^ f 
\ Lm 0 


si X dm 


rrr 


+ d. 
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m 


/ 

Lm„ 


dm 

— T 

nr 


r m 


+ d. 


11 


u m 0 


f dm 

J m 


In order to evaluate the integrals involving si X in the above expressions, 
four cases must be considered. These cases are at the top of the next page. 
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1 . = 0 but dg > 0 

2 . d 7 > 0 but dg = 0 

3. d 7 > 0 and dg > 0 

4. <17 = 0 and dg = 0 . 

It should be noted that d 7 and d 9 can never be less than zero because 
d 7 = [A 0 ] T j^u 0 + u 0 | 

and 

, -1 rA ,Ti 

dg = ^2 l A (P u 0 

Also, if either d 7 or d 9 is zero, then d 8 is zero because 



and if d 7 or d 9 is zero, then cTj is zero or dg is zero, respectively. 

For the first case, since d 7 = 0 and d 9 > 0, the expressions for " 7 ( 10 ) 
and y(m.) become 



The integrals appearing on the right -hard sides of the above expressions are 
now easily evaluated to give: 


il 



T(m) = v 0 - i |(^- + b ») f "(^) + So (m - m « ) } ' 

y(m) = Yo - S {(^ + d..)'"(^-) + d.o(^ - £) j • 

To complete the closed-form solution for this first case, it should be 
noted that 


dx _ dx / dt \ 
dm dt 


1 

m 


(©■(-i)S-(-i)' 


v 0 - 

m 


dx 


- m 0 ) 


This expression for can be integrated to give 


i m 

js(m) = So -(^■)' r o / dm 

mo 


(®)[( 


, , V m 

nTt B + Eoj / (fa m - in mo )dm 
^ ”• 'm 0 



Thus, the closed-form solutions that approximate the entire system of 
ordinary nonlinear differential equations when dj = 0 and dg > 0 can be cal- 
culated for any t > t 0 as follows: 
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m(t) = m p - m(t - t 9 ) 


l(t) ~ + Xq (t - to) , 

ri(t) = n 0 + tr 0 (t ~ t 0 ) , 


v(t) ” To - 


v(t) = ?o - 


TO 


TO 




+ ^°) in {m^) + to 0* - m 0 ) 



' 


*(t) = X 0 + ? 0 (t - to) + g 0 

+ (^)(^f + 5,)[“ (~) - (» - m o ) 


For case 2, since d T > 0 and dg ~ 0, the expressions for v{m) and 
y{m) become 


T,n) -* -sl^r 


m 

/ 

L m » 


dm 

m 


+ bn 


m 

f 

mg 


dm 

m 


So 


m 


/ d ™ 

m 0 


T(m)= Y 0 


m 

/ 


dm 


+ d 


10 


TO 


f 

Mr 


dm 


+ d 


li 


• m 

I 

M n 


dm 

m 


The integrals appearing on the right-hand sides of the above expressions 
also easily evaluated, and the results follow: 


are 


V(m) = Vo - 4- 

m 



+ ^°) (m^) + " m °^ 


y(m) * yo - -r 
m 


+ d 10 )f“~- - — ) + dji in 

' m o to / - \mo J 
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Then, to complete the closed -form solutions for this case, it should be noted 
that 


dx _ d>T / dt \ _ / 1 \ dx _ / 1 \ _ 
dm dt \dm ) ~\xn/ dt \ rn/ V 



d3? 

This expression for - — can be integrated to give 


■Jt(m) 



+ (-^2) + / (in m - in m 0 )dm 

/ L V 7 'mo 

m 

+ go / (m - mo) dm J 


m n 


= X° - (ij v 0 (rn - m„) 


b ")[ mi "(m 0 ) ' ( 


m - m 0 ) 


Thus, the closed-form solutions that approximate the entire system of 
ordinary nonlinear differential equations when d 7 > 0 and dg = 0 can be calcu 
lated for any t > to as follows: 
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m(t) = m 0 - m(t - to) , 

\(t) = A-o A-o (t — to ) , 

u(t) = G 0 + t 0 (t - to) , 


v(t) = Vo - + <W(^ - ~) + d u «n (^) ] , 

v(t) * T, - i[(^ + bo) (—) + lo(m - “o)] . 

!i(t) = s 0 + 7 o (t - t 0 ) + go [ — ~ 2 "^ ] 

+ (sOfalr + b °) [ m in (37) - (m - m » > ] • 


For case 3, where d 7 > 0 and d 9 > 0, the expressions for v(m) and 
•y(m) remain the same as the ones given on page 11; i.e., 


v(m) = v 0 


1 /- 

rm . 1 

r dm 

m 

+ f 

dm - 

r rn 1 

r dm 

m dl 

J 37T 

n + b » 

m 

l 

Lm 0 J 

m 0 


Lm 0 J 


+ So 




t 


V (m) - To - i I 


m 


m r 


dm 


" m 
r 

\/ X dm 


’ p dm 

J 

— — Y — 

m' 

+ d io 

J sr 

-m 0 



L m 0 

m 

r 

dm 

m 
, r 

Vx dm 

. J 

nTX 1 

and J 


mo 


m o 



+ d 


li 


m 


i-m 0 


r dm 
J m 


expressions are evaluated with a table of integrals. This gives 
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To facilitate writing the above expressions, let 
<f> = 2rTd 7 nTx + dgm + 2d 7 
£2 = 2\Tdg + 2dgm + dg , 

and $ 0 and £2 0 will denote $ and £2 evaluated at the initial conditions. Then 
the expressions for 7(m) and y(m) can be written as follows: 

V(m)= v„ . [‘n (i) - in (;=-)]. 

+ ^(^ar) <n (£;) + 5 » Jn (-^y *• (m - m « > } ’ 

+ ^ (W") + d >»(^- - ST ) + d » ,n (^) } ' 

To complete the closed-form solutions for case 3, it should be noted 
again that 

dx _ dx / dt \ _ / 1 \ dx _ /J_\ _ 
dm dt \dm / \m/ dt \m/ 
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Thus, 


dX 

dm 


= -(sH + (s>)| a *(-7T)[ in (£) *'"(^)] 

+ *(*%) tn (It) + 5 ° (57) + ? * (m - ra °>} 


dx 

This expression for - — can be integrated to give 

dm 


x(m)= x 0 - (i) v 0 / dm + 4j | - 
' ' m n ' 



m 


(in $ - in $ 0 ) dm 


m» 


m -i - m 

-/ m - in m 0 )dm + -A- J (in SI - in S2 0 )dm 
m 0 -1 ^ m 0 


m 


m 


+ bo / (in m - in m 0 )dm + go / (m - m 0 )dm 
m 0 m 0 


The only difficult integrals to evaluate in the above expression for x(m) 
are the integrals involving in $ and in fl . The following steps will explain the 
m m 

integrals f in $ dm and f in dm. Both integrals will be evaluated using 
mo mo 

the integration-by -parts formula; i. e., 

f u dv . = uv - f v du 


m 

For f in $ dm, let u = in $ and dv = dm. Then, 
m 0 


m m 

f in $ dm = m in $ - m 0 in $ 0 - f 
mo mg 
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Then note that 


_ / 1 \ / 1 \[" \/~ck do \f X + 2d 7 d 8 m + d 7 do 

= w fen?) [ — -rsr* — 

•• 



^8 1 + f ^"^7 \ 2 ^7 

2 *T5; Jx v^/ ‘^"rndgNTx 


j_ 

m 


n/ d 7 

nW"x 


When this result is 
the result is 


m 

substituted into the expression for J in $ dm, 

m 0 
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m 


m 


f in $ dm = m in $ - m 0 in $ 0 - J m 
mo mo 

m 


J_ nT d* 

m 


dm 


m 


dm 

sTx 


= m in i - m 0 in $ 0 - J dm + nTc J 

m 0 mo 

«. m in * - m« In * 0 - ( m - m 0 ) + ^ (~j 

m 

A similar procedure is used to evaluate f in A dm; i. e. , 

mo 


Then 


m 

f in S2 dm = m in J2 - m 0 in O 0 

m„ 



_ /A*nT"[' 2s ^fl + 2dgm + dg 

" ^ L m j 

NfdT 

— *■ n 

nTx * 


m 

Now, as before, this result is substituted into the expression for J in S2 dm 
to give m 0 


m 

m 1 

r K r 

hTI 

J in 0 dm = m in 0 - mo in £ 2 0 - J m 

M 1 

TT 

X 

m 0 

mo 1 




dm 


m 


= m in £2 - mo in £ 2 q - ^"dg f 


mo 


mdm 

7TT 


19 



m 


For the expression f » a table of integrals is again used to give 


mo 


r 1 mdm _ \f~X \Tx n 
4 <*• " d*) 

m o 


- (&) “ & 
N 7 m n 


= - (^)ln (i)] 


Thus, 


J in Q dm = m in Q - m 0 in n 0 - (t?) p* - ^ 

.-(A 8 ) in (t)] • 

With these results, the evaluation of the integrals appearing in the expression 
for x(m) can now be completed; i. e. , for case 3, where d 7 > 0 and d 9 > 0, 

X(m) = %, - - m,) + f 0 < m J 

+ (4,) b t [min - (m -m,)] 

+ .(s)(^D ( m [ ,n (£) - ia fe)] 

) iu (^) } + (a?) { m (£)■ 
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Now, again, the closed-form solutions that approximate the entire 
system of ordinary nonlinear differential equations when d 7 > 0 and d B > 0 can 
be calculated for any t > t 0 as follows: 

m(t) = m 9 - m(t - ) , 

\(t) = \ + %(t - 1*) , 

u(t) = Uq + ir 0 (t - * 0 ) * 

? (t > - 1 j-(f -t?) - (fa) [ ta (£) - In (^)] 

+ ^i O (^) + d 10 (J 7 -i) + d 11 in(^)[ . 

v(t) + *> <‘ -«.) + &■ { -fa) [ in (if) - *> (^)] 

+ (^) in (^) + s » fn (^-)j - 

x(t) = Xg + V 0 (t - t 0 ) + g 0 + ^4^ b 0 |m in 

- (m - m °>] -(&) (fcn?) H® - in (5)] 

+ fa ) ® } + (I s ) { fa ) in (£) 

-(i)[fa-^-(fa)‘°(%)]} ■ 

For case 4, when d 7 = 0 and d 9 = 0, the aerodynamic forces and thrust 
forces are zero so that motion is two-body motion for which closed-form 
solutions are developed in Reference 2. 
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The preceding work completes the closed-form solutions for the ascent 
case which can be used to approximate the optimal atmospheric motion of a 
rocket-powered space vehicle. Some of the results will be used in the next 
section, which is concerned with the closed-form solutions for reentry motion. 

Reentry Case 

For the reentry case, F is assumed to be zero and m is assumed to 
be a constant equal to m 0 . Also, for the reentry case, q, which is a linear 
combination of the instantaneous stagnation point heating rate and deceleration, 
is defined by the following differential equation: 



and the initial conditions qo = 0. To simplify the writing of the equation for 
q, define 



and 



+ 4q c c 2 

A Xj 

a 



Then 


q = 


h + 


h T 


p 

I pi * 


Since F = 0 and m is a constant, the equations of motion for the reentry 
case become 
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x= V 


9 


- [A] JL 
Ipl 


+ b + g 


9 


where 



g = 


GMx 

R 3 


( 


Now, as in the ascent case, the multipliers \(t), u(t), and y(t) (where 
^ s t ^ t ) are introduced to allow the following expression for H to be 
written: 

H = "X T u +u T |[A]^ + b + gj + y ^h . + h T y|0 . 

As in the ascent case. 


4 T 9 H IT 9 H . . 9 H 

* = ~bT ’ u - -aT- andy= -R ' 

Now the quantity J, which is to be maximized or minimized by the 
choice of p(t) (where ^ ^ t ^ t ), is defined. As in Reference 1, the choice 

j = is made, and it is desired that J be minimized by the choice of a 

control function p(t). Thus, p(t) must maximize H at every t such that 
t 0 < t ^ t^. An examination of the above expression for H will indicate that 
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T - 

p = [A] tT + y h 


is sufficient to maximize H at every t, where ^ t ^ t^. 

When the above result for p is substituted into both the equations of 
motions and the multiplier differential equations, the entire system of simul- 
taneous ordinary nonlinear differential equations for the reentry problem is 
obtained; i. e. , 


x = v 


v = [A] 


[ A] ^u + yh 


7 (u T [A] + yh T ) ([A] T n + yh) 


+ T> + g , 


— T 

q = h + h 


i T = 


[ A] T u + y h 


7 (u T [A] + y h T ) ( [A] T 1+yh) 


T -T 

a [A] + y h 


7(u T [A] + y h T )([A] T u + y H) 

T / 9 B a 9 g\ (dh\ 

\9x 9x/ ~ y \9Y/ ’ 


adAl V) x v d ± 


9x 


9x 


_T 

u = - 


u T [A] + y fi T | ( 9([ A] T n) 9h_' 

| 9v 9 v 


i7(u T [A] + yTi T ) (tA] T u + y h) 
_T^9b^ ^9h^ 


y = 0 


This system of ordinary nonlinear differential equations can be solved simul- 
taneously by numerical integration if the initial conditions Xq, v 0 , X 0 , Ho, 

and y 0 are given; or, as in the ascent case, the quantities, X, u, h, h, [A], 
b, and g can be expanded in truncated Taylor Series about the initial conditions 
which will allow numerical integration by quadratures to be used. 
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To continue and obtain closed-form solutions for the reentry system of 
differential equations, the following assumptions are made: 

^■(t) = + ^o(t - to) * 

n(t) = Ho + %(t - to) , 

h(t) = "Ho + "Mt - to) , 

t A (t ) ] = [A 0 ] , 

• • 

b(t) = b 0 + 5 0 (t - to) + &-(t - t 0 ) 2 + ... 

• » 

g(t) = go + g 0 (t - to) + ^L(t - t 0 ) 2 + ... 

• • 

h = ho + h 0 (t - to) + . to) 2 + ... . 

As in the ascent case, any number of terms can be used in the Taylor 
Series expansions for b(t), g(t), and h(t); but only the constant terms will 
be used in illustrating this approach. With the preceding assumptions, the 
equations for v and q become: 


lA 0 ] ‘in, + ^(t - t*)J + ylfij + E„(t - to)J 


| -i([o« T + - *<>)] CA 0 ] + y[t7 + t, T (t - t,)]) ([A 0 ] T fa + Ub(t - t„)] + yfE, + h„(t - to)]) } 

T • T 1 

h 0 + K 0 (t - t 0 )J {same term as above in the v equation} . 
To simplify the expressions for v and q, define 


d t = [ A 0 ] |[A 0 ] T [u 0 - Uot 0 ] + y(E 0 - iotoll , 
d, = [A 0 ] | lA 0 ] T Go + yEq } , 

d 7 ={[u 0 - Uot 0 ][A 0 ] + y[hg - h 0 t 0 ]^ ^ [ A 0 ] [HJj - ^ot 0 ] 
+ ylE 0 - Kot 0 ]| , 
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= 2 - n^t 0 ][A 0 ] + ylH^ - hjtojj { [A 0 ] T u 0 + yh 0 } 


d s = {u^IAq] + yhj^j | [ A 0 ] "tig + yh 0 


} • 


-T -T ( 
d io = l h o - h 0 t 0 ] ^ 


T _ j_ — - \ 

[A 0 ] [uo - Ugtg] + y[h 0 - h 0 t 0 ]f 
r * ,T. , .A . fTf r A ,T r= • 


d n ~ [fro - hg t 0 ] ^ [A 0 3 u 0 + yh 0 | + h 0 [A 0 ] [ug - Ugt 0 ] 
+ y[h 0 - H 0 t 0 ] | 


= Ho T { 


T * * 

[A 0 ] + yh„ 


} • 


Then the equations for v and q become 


d, + dot 


d 7 + dgt + dgt 2 


+ + So 


q = h 0 + 


dm + d tl t + d^t 2 
\l d 7 + dgt + dgt^ 


These equations can be integrated in closed-form with the aid of a 
table of integrals, but four cases must again be considered as in the ascent 
case. These cases are 


1. d 7 = 0 but dg > 0 

2. d 7 > 0 but dg = 0 

3. d 7 > 0 and dg > 0 

4. d 7 = 0 and d 9 = 0. 

Also, as in the ascent case, d 7 ^ 0 and d 9 ^ 0; d 7 = 0 implies that both 
= 0 and dg = 0; and d 9 = 0 implies that both = 0 and dg = 0. 

Thus, for the first case, 

i 3, r _ 
v = -J- + b 0 + g 0 

Vd 9 
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* _ Vi + d 10 + + d 12^ 

q - n 0 f _ 


Nfdlt 


These differential equations are easily integrated to give 

v(t) = V 0 + (— i- + B 0 + (t - to) 




q(t) = q 0 + 


(h, + iSL.\( t -t 0 ) + In (l) + ^l(LzA\ . 

\ W ^ W nT*' 2 ' 


Then, since If = v, 


-^2- + b„ 
\Td 9 


x = v 0 + j — — + b 0 + g 0 (t - t 0 ) 


which is easily integrated to yield 

x(t) = X 0 + v 0 (t - t 0 ) + + b 0 + go ^ . 

Thus, for case 1, the entire closed-form solution which approximates 
the original system of ordinary nonlinear differential equations is given by 

Mt) = + ^o(f - ^0) > 

u(t) = u 0 + u 0 (t - t 0 ) , 

Y(t) = 7 0 » 

v(t) = v 0 +(A_ + b 0 + g 0 ) (t - to) , 

J 

x(t) = X 0 + v 0 (t - t 0 ) + ( -4--+ b 0 + g 0 ^ 

^d" 9 


q(t) = q 0 + 


N - g + i« fn (f\ + in • 

, -TS,/ -fi, w 2 > 
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For case 2, since d 7 > 0 and d 9 = 0 , 


V = — i- + b 0 + g 0 

n/~c£ 

q= h 0 + d iQ + d U- t-lu.?. . 

These differential equations are also easily integrated to give 

v(t) = T 0 .+ ( -A. + b„ + g 0 Vt - t,) 

j 

Then, since x = v , 

It = v 0 + [_A_ + b 0 + g 0 ) (t - t,) 

y-rs; ) 

which is easily integrated to give 

x(t) = xi + ? 0 (t - t„) + (-A- + b, + go ' — — 

Then, for case 2, the entire closed-form solution which approximates 
the original system of ordinary nonlinear differential equations is given by 

X(t) = \ + \(t - to) , 

o 

u(t) = u 0 + Ti 0 (t - t 0 ) , 

r(t) = r 0 , 


v(t) = + [ — - + b o + go) (t - to) * 

\Td 7 
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*(t) = Xo + v 0 (t - to) + b 0 + g 0 \ (t ~3 )2 , 

VnTX 


q(t) - 4. - ( h, + - i,) * + ^ (~^) 

\ -J *7 var v 2 ' 3 ; 

For case 3, since d 7 > 0 and d 9 > 0, the original form of the 
differential equations for V and q are unchanged; i. e., 


d t + a^t 


\/ d 7 + dgt + dgt? 


+ T>o + g 0 


q = hp + d l£ 1 d 12^ 

^1 dj + dgt + dgt? 

As was done in the ascent case, define X = d 7 + dgt + dgt 2 . 

Then the differential equations for V and q can be integrated at least symbol- 
ically to give 

?(t) = 7 0 + j, /-£-+ a* (B„ + fo )(t - 1 0 ) 

t(, N fX t 0 -v/X 

q(t) = q 0 + h 0 (t - t 0 ) + d I0 + d tl /-^+ d^ f-^~ . 

to \ fx t 0 *Tx to *Tx 

With a table of integrals, the expressions involving *J~X can be integrated 
to yield 


f _dt_ = J_ in / 2^ nTx + 2dpt + dg \ 
t 0 Nfx nT^ V^/dg + 2dot 0 + dg/ 


r tdt _ _1 

l 7 z~ di 


•Jx.-Tx, -(js-\ln( jjjLLA 

2^ dg J y 2nZ dg nTX o + 2dgtg + dg ^ 


29 




Let S2 = 2sfd^ 'Vic” + 2d g t + dg and S2 0 be the value of SI at t 0 . Then 
the expressions for V(t) and q(t) can be written as 



Since x = v, 


*(t) = r 0 + V 0 (t - to) + -A- f (In St - in S2 0 )dt + (b 0 + g 0 ) ^~^- 
+ 4 ) f (m fx _ \Tx 0 )dt - — f (in 0 - in 

* \i 2nTH 7 t 0 

m 

The integral J in S2 dm was evaluated for the ascent case; and when 
m o t 

m and mg are replaced by t and t 0 , the integral J in ft dt can be evaluated. 
That is, 4 
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J in fll dt = t in - t 0 in f2 0 - 

to W*/ 


Also, using a table of integrals, 


*nc - ■m: . ( - 


in 


\2nTX| 


© 


( sfx dt = ( 2d 9 t + (2dgto + de)^ / 4d^ - _df\ r 

J tj 4de 4dg \ 8de ) J 


dt 

\Tx 


= WIT- trfvTXfl + dfi (\T3T _ */3Q 

2 4dg 


Thus, 


x(t) = + V 0 (t - t 0 ) 


d« ( . / SI \ (l 

+ _ _ L. J t in ( — ) - / — 

1 v 0 "/, 


sTar/L 


nTx - \Tx^ 


- / AJU (AjU + A_ . ^x {t . t ) 

v n »/J 3 d . ( 2 

* ( ^)(^)» ( a 

(£)]} ^ • 


jL-i n 

2^ 


Then, for case 3, the entire closed-form solution which approximates 
the original system of ordinary nonlinear differential equations is given by 

^■(t) = + N)^t - to) , 

TT(t) = Ho + cr 0 (t - to) , 
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y(t) = To 


f 


r(t) = v 0 + A_ln (§) + & 

■r^ vv d » 

+ 05o + lb) (t - to) , 

x(t) * x'o + V 0 (t - to) + 


-fx Jjx; . ( A_\ <n 

W^/ 


(4) 




in;) 




+ A 

d Q 




+ / 4d 7 di) - 3dfl 


8d 0 


feH4) -(*)"■ (4)| 


q(t) * q 0 + h 0 (t - t 0 ) + A_f n 


\Td^ 


. Sj 




A 

2nTc^ 

-i? 


in 


(§)] + 4*p-v/*-- 


+ ( 3 . d 8.^. 4 .^ j in|~ 
\ 4d 9 '/c^" / 


f2 n 


For case 4, when d 7 = 0 and dg = 0, the aerodynamic forces are zero so 
that again the motion is two-body motion for which closed-form solutions are 
developed in Reference 2. 


The preceding work completes the closed-form solutions for the reentry 
case. These closed-form solutions can now be used to approximate the reentry 
motion of a space vehicle. In the next section, the reentry and the ascent 
closed-form solutions will be used to develop expressions for the closed-form 
partial derivative transition matrices. 
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THE PARTIAL DERIVATIVE TRANSITION MATRICES 


In order to define most easily the partial derivative transition matrices 
for the three ascent cases and the three reentry cases, a new vector of dimen- 
sion 14 is defined; i. e.. 



Then, for each of the three ascent and reentry cases, the partial derivative 


transition matrix is the 14 by 14 matrix 


"BZ 

“ L 8 Z o. 


This matrix can be eval- 


uated at any t, where t 0 =£ t ^ t^, by explicit differentiation of the expres- 
sions for x, v, q, X, u, and y with respect to xj>, v 0 , q, T,,, uj, and y 0 

9 Z 

Also of interest for later use is the vector -t—. The transition matrix 

at 0 

9 Z 

and the vector ~ can be written more explicitly as follows: 

dto 


' az ' 

ix — — 

az 0 


ax 

dx 

ax 

ax 

ax 

ax 

9*0 

av 0 

aq 0 

aXo 

au 0 

ayo 

dv 

8v 

dv 

3v 

dv 

dv 

8X 0 

a? 0 

aq 0 

ax 0 

0U O 

ayo 

JS L 

JL 

3q 

ag 

aq 

aq 

ax 0 

dv 0 

. aq 0 

ax 0 

9Ug 

ay 0 

ax 

ax 

ax 

ax 

ax 

ax 

dx 0 

av 0 

9q 0 

axo 

9u 0 

ay 0 

au 

3u 

8u 

85 

3u 

an 

ax 0 

av 0 

aqo 

ax 0 

0Ug 

ayo 

J2L 


j*L. 

JhL 

iJL 


_dx 0 

dv 0 

aqo 

ax 0 

au 0 

ay 0 
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dx 
9 to 

3T 

a t 0 


az 

9 to 


_9q 

at 0 

ax 
9 1 0 


9u 
9 to 


9y 

9 to 


' y ^ 

To illustrate how the different components of -r-=- amd — - are obtained, 

dZ 0 oto 

the expressions for ■ and will be derived in detail for case 3 of ascent 

9x o 9 to 

flight. The other components (including the other cases for ascent and reentry) 
can be derived in a similar manner, but their derivation will not be shown in 

detail. To start the derivation of -r=- and •*— , note that for case 3 of ascent 

axo ato 
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To ease the problem of writing the above expression for x, define the scalars: 

c * ' [ m ,n (t) - (m - m °’] G? ? ) • 




i2\ r dly 


Then 


and 


x = x 0 + v 0 (t - t 0 ) + fo [ ^ t "~2" i) ' ] + c l b 0 + c 2^l + ^ 

II = - (f 0 ) - (H) (It) ^ tf) 

+ 5 *(it) + c ’(lf) +a * (l|) • 


To complete the derivation of — , the unevaluated terms in the above 

OXq 

expression are taken one at a time. Since 


_ _ -GM 
g ° .T Z \% 


^7. 


(xq 5^0 ) 


then 


GM 


-GM 3GM Xq x. _ 

9x ° (xjf Xq) 2 (X<F Xq) 2 (55? ^o) 2 


T 

3Xfl Xq _ 1 

.(xf *b) 


Similarly, 


b ° = -(-T 1 ) l7 r' (°A + 
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where p, v , c , 17, and c are all evaluated at t$ and thus depend on "x^. 
Then, r A a 

^ -(t) 'V ( C A + ^ L J) \(^) 

- (pA r)( C A + 2,,C Lj)( T r 7 r) 9^ 




|7 | /c. + 2ijc_ 
r l A L 


av /pA \ 

J) 3^ “ \— ) lv r' V r 


dc , 


^ ,2 


+ 2 -r— C‘ + 2r]C -r=^ 

1 9%, L L 9x 0 

a a 


Since evaluated at t Q is given by 


d .2i 

dx n 


v r = v 0 - to X Xo , 


then 


9 v 
i 

9x n 


0 - co 3 a; 2 


00 3 0 - GO j 


—CO 


2 0 


where co t , co 2 , and co 3 are the components of the constant vector co explained 
in Reference 1. Also explained in Reference 1 are subroutines for the atmos- 
pheric functions and the aerodynamic functions. The subroutine for the atmos- 
pheric functions yields values for the density p, the pressure P, and the 

velocity of sound v as functions of altitude and also the derivatives -~r , 

S dalt 

dP dv s 

, and • The altitude at is given by 


dalt 



where is assumed (for this development) to have a constant value equal to 

the mean value of the earth's radius. The subroutine for the aerodynamic 
functions yields values for the coefficients c. , c , and tj as functions of 

A 1j 


a 


9c, 


3c, 


a. 


Mach number M and the derivatives — 

3M 3M 

Mach number M is computed by the following formula : 


, and 


J2L. 

3M 


The 


M = 


9c t 

dp dC A 3 ^ a 

Now the quantities -r=~ , ■ , — - , and —r — — appearing in the expression 


9xo 3x 0 31J0 


for -—I 2 - can be written in more detail, but first the quantites \ ^ and 

3xo 9x 0 


3xo 

must be written. This gives 




0 

This completes the derivation of . 
— — ) Cj does not depend on Xq so that 

9"Xrt 


Now note that (in the expression for 

9 C 

-£=*• = [0, 0, 0] . Then to continue. 


dt = [A 0 3 [A 0 ] T [uo - . 

Thus, 



To evaluate the quantity 



a[A„] 

3x5 


[u, - Uo(^f) 


note that 


Then, 


- *•&)]} = [* - Mt)J (r) 1 V \ (U) + !p V \ v r {S + (V) ‘V 3 -sir 

■ !'![“. - '©)]} I’, (f ) \ (■ - »x) ft) ‘ (£)(■ - -x) - w r ) \(#> 


8F 

The only quantity not already derived in the above expression is -r=- . From 
Reference 1, ^ 


F = F g + A e (P 0 - P) 
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where F , A , and P 0 are constants and P is obtained from the atmospheric 

o G 

subroutine. Thus, 


9F / SP_\ ( 9ait\ 

9x 0 " e \ Salt/ Sxo/ 


This completes the derivation of the expression 



To complete the expression for , the expression -1^ 

OXq OXq 

evaluated. To do this, it should be noted that 


must be 


_lT -T 1 

U = -A - 

u t [aj a([A] T u) t 

f— ) ! 

m 

|tt t [a] 1 av 

w j 


Then 



Note that 


and 
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Now the above expression, along with the initial conditions, can be used to 

evaluate TT 0 , and then can be computed; but before attempting this task, 

3Xo 

the matrix [ B] and the scalar c 4 will be defined as follows; 





Then, 


Uo = -Xo — E B 0 ] 


[A n ] T Uo 
I [ A 0 -] T UqI 


+ 



v _ + c 8 0 "0 
r 0 
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and 


aun 




a[A 0 l ) 

' u ol 


ax n 


9x n 


|[A 0 ] T u 0 


l[A 0 ] , u 0 


-ft) 


(\*) 


v ^ U 

r 0 \8x 0 / 


C 4 n 


(V°) + (\ "«) 


'9v 

v ax 0 


{ d [ B ] ( lT 
■ d _ 11 ([A 0 ] u 0 

dx 0 \ 

j — ) u n / , , and have not yet been calculated. Before 

(\ 9x o / °) ’ £E 0 8c 0 

attempting to calculate these quantities, some more constants will be defined; 
i.e., 


■& 


f 9c, 


2c + M 
L 

a; 


or 


8M 


° s = Isr/ c l 


1 - 2tjc. 


a 


a , 




• (&) 


'9c, 


1 - 4tjc 


a 


a 


3M 


- 2c 2 


(— ) 

'a W 


Then 


[B 0 i = c 5 (v r u 0 T ) - c £ (vj' u 0 j I - 


( V r U °) 


v r v s 


( v r v r) • 


where c 5 , c 6 , c 7 , v , and v are all evaluated at t 0 .and thus depend on x 0 . 
r o 


Now the quantity 


ity | (lA 0 ] T u 0 ^ j can be written; i.e., 
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To complete the evaluation of the above expression, 
be written. This gives 





must 
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Note that~rr ?2 and — 777? ' are assumed to be available from the aerodynamic 


subroutine of Reference 1. To continue evaluating 
will now be evaluated. Since 


{SSt) 

\3V’ 


the expression 



then 
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8u f 


To complete the evaluation of 7— I 

3x n 


f 



45 



3 2c a 8c, 

Again note that needed to compute is assumed to be available from 

the aerodynamic subroutine of Reference 1. This completes the evaluation of 

-rr“ and thus also ■£■=*• • To continue the evaluation of ,, the next term to 
dxQ dxo dx 0 

be evaluated is 4=?“ . Since 
axo 


then 


°2 = 


1 \ m 

-T 9 




S2n 


»a. -4„ 

S’Xo m 2 


2(n T d^) 


+ J- hi) 

V5 ;\ a ) 


2— fi„ (*) -J^\ 

2(-Ti,) 3 \ 8x V L v V \ m °/ . 


) '"(ft)! 

ill 


/ ad 

7 l 8^ 


$ 


•(*) 


axo 


8^1 




a 


d 7 * [ ^ 0 T + ^ T (^f)] d l 


and thus 


ad 7 

3xn 


_ T : 
Uq + 


^)]i + (^K(iS 


To complete the evaluation of 4^ » the terms 4=?- , 4=^ » • ? , and 

3XJ, 3x 0 axQ ax 0 

(|-) 

-z-S — must be evaluated. First, 

OXq 
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Then, since 


d 9 = + ^ » 

fie (*)[*»)• <5)1 • 

To evaluate 4^2- , note that 
9 xo 

d^= -[A 0 l[A 0 ] T n 0 (i) . 


Thus, 


Note that the expression for j — ^u 0 j 

( T ) 

for \ 8[ >L u 0 ( with the exception that u 0 is replaced by u 0 . This 
( dl ^o ) 

completes the expressions for 4^- and -r— ?• . Also, since 

dX() ox 0 


is the same as the expression 


■ -(*> 


• T - 
% d l , 


then 


•ft) *fe) 

which will be needed in the expressions for y ZS-‘ m and — ■ “ ■■■ ? ■ 

cJXq dXfl 

$ = 2nTc^ \/"x + dgm + 2d 7 


Since 


and 


= 2\Td^ M Tx + 2d 9 m + dg , 
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then 



\Qn) _ / 1\/9Q\ Q /afl n \ 

r . , 


S2 0 


[■/rW [**>) \ 


£&' 

8x 0, 


wi(^)*t(^)*-(«)*(3 


3X 3X 

In the preceding expressions, the terms ■%=■ amd -=?• are evaluated as 
follows: dX ° dX ° 


S ■ H * -(>!) * "'(Si) 

-(a) ■ 


This completes the evaluation of the term The next term to 

OXq 

evaluated in the expression for is 4=^ which was evaluated already 

8x5 3 X, 



9 Co 

order to obtain . Thus, the only term left unevaluated in the expression 

dXfl 


for -7- is the term 4=2- . To evaluate 4^?- , note that 
9x 0 9 xq 9x 0 


c 3 = 


/ 1 \ m . / n \ / 1 

i 3 '/ ^ \ a ») ' l* 


- (> /x - "/x 0 ) - 


2^/ \ 


iL 


Then 


|2- = (4*) 

Xq \mV 


£ 

3X0 


m 


/'£i\ in / iL\ + /J2_\/'£A| 


2(nT^> 


_J_ /MAI ..fa/O-Nf— (^\ ^ / ad » 


9xo 


( _A_\ K) 

K) 1 


ax o 


> This completes the detailed derivation of . A similar approach is 

dXo 

Ox 

used to derive an expression for — . That is, since 

3 to 


x = Xo + v 0 (t - to) + g 0 


(t - t 0 ) 2 
2 


+ c t b 0 + + c 3 d 2 


then 


9x _ _ 

3tn " V ° 


S» (t - *»> + + (If) ^ + ( 1^) 53 


0 b 0 d 

Note that — l , and -^f 2 - are zero because b 0 , d lt and d 2 are constants. 

9 to 9 t 0 9 to 
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The terms 




and 



must now be evaluated. 


Since 



then 



Note that m = m 0 - m(t - t$) so that = m. 


Thus, 



Next, since 




then 



Note that 

$ = 2Nfd^ + dgm + 2d. 7 , 

£2 = 2hJ~d^ -sfx + 2 dgm + dg , 


X = d 7 + dgm + dgm 2 
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Thus, 



This completes the expression for *rr* . 

dtg 

The preceding work has shown in detail how to compute —=~ and 


The other terms in the partial derivative transition matrix -r— — and the vector 
dZ 

-rr- can be written in a similar manner and, in fact, most of these terms are 
dto 

much easier to write than the ones used for illustration purposes. In the next 

0 Z 

section, the partial derivative transition matrix — — will be used in deriving 
a variation-of -parameters integration scheme. 0 


VARIATION OF PARAMETERS INTEGRATION 

For more information about this approach to variation of parameters 

Integration, see References 2, 3, and 4. Variation of parameters integration 

is to be used to obtain Z (t) where Z (t) is the vector 

a a 
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z (t) 
a 


x(t) 

a 


q 0 ! (t| 

a 


X i(t 

a 


*«« 


defined by the differential equations as given in Reference 1 (also listed in 

a previous section, Development of Closed-Form Solutions) and the initial 

conditions Z (to) = Z 0 . 

& 


Now define Z = Z - Z, where Z is the approximate vector 

pa 


differential equation used to derive either the ascent or the reentry closed- 

form solution and Z is the actual vector differential equation mentioned 

a 

above. For example, in the ascent case 


v 

(!) + (IK + io 

I [ A 0 ] [v^ + ^(t - t 0 )] ! q 

-m 

u 0 

\ l[A 0 ] T [u 0 + R 0 (t - to)] I + [uo T +t^(t - to)]b 0 
4 
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and 



Then Z = Z + Z . Now a vector Z’(t) is defined to be a set of values 
a p 

which, when used as initial conditions in the closed-form solution, Z(t) will 
make Z ( t) = Z ( t) . The following diagram might help to understand the 

3 . 

definition of Z'(t) . 


lt,Z(t) - Z,(t)J 


The diagram illustrates that Z (t) is obtained by numerical integration of 

a 

the actual differential equations with the initial conditions (to, Z 0 ) . Also, 

a value Z(t) = -Z Ct) is obtained using the closed-form solutions and the 
a 

initial conditions (fo, Z ' ( t) ) because this is the definition of Z '( t) . Given a 
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value of Z (t) a value for Z'(t) is obtained from the closed-form solutions 
a 

by letting Z (t) = Z 0 and interchanging the role of t and t<,. To be more 

3 

specific about the definition of Z’(t) let the closed form solution Z(t) be 
denoted more generally as 

Z(t) = g(t, to, Z,,) . 

Then from the definition of Z*(t) 

Z (t) = git, to, Z * < t) ] 

3 

or 

Z'(t) = gt to, t, Z (t)] 

3 

since the closed-form solution works backwards or forwards when t is 
replaced by to. From the expression 

Z'(t) = glto, t, Z (t)J , 

a differential equation for Z'(t) can be written. This differential equation is 
the variation-of-parameters differential equation, and from the diagram the 
initial conditions for the variation of parameters differential equation can be 
seen to be Z'(to) = Z 0 because Z 0 used in the closed-form solution over a 

zero length of time will still be Z 0 , which is also Z (t 0 ). To write out the 

d 

differential equation for Z ’ (t) explicitly, note that 



Z = Z + z 
a P 
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Thus, 


a(* ■<*>)■ 


Sgltg, t, Z^(t)] j &g [to» t» (t) ] 


8to 


at 


f aglto, t, z a (t)] | . 

4 . 

a g [to, t, z (t)] 

cL 

j az”Ttl j z 


az (t) 

' a / 


a 


Note that Z (t) is used to evaluate Z and Z . Now it will be shown that the 
term a P 


9 g(t 0 , t, z (t)) a g [to, t, z (t)] 

a + a + 1 

9 gft 0 , t, Z (t) 1 ) . 

ato at 

az (t) 
a 


= 0 


so that 


Tt - 


a g [to, t, z a (t)j j . 
az “(t) [ ^ 

a i 


To see that the term just mentioned is zero, note that the initial conditions for 
the closed-form solution are constants. That is, 

Z 0 = g[t 0 , t, Z(t)] 

for all [t, Z (t) ] on a particular closed-form trajectory. Thus, 
d Z f 


dt 


= 0 , 


and using the above expression 


dZ^ 

dt 


■ 8g(to, t, Z (t) ) + agltp, t, Z (t) 1 

at„ at 


9 gttfl, t, z (t) ] ^ 

— sz(i) ' z 
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Now if a particular Z (t) is considered to define a closed-form solution, then 
Z (t) substituted in the right-hand side of the above expression will give the 
result that 


8g(to, t, Z a (t)) 8g[to, t, Z a (t)] 


9% 


at 


9g[to, t, Z a (t)] J # 

az (t) 

a 


= 0 


Thus, the variation-of-parameters differential equation can be written as 


I < 2 '<*» 


3g[t,, t, Z (t)J | . 

- — — | z p • 


Now note that since 

Z (t) = g[t, to, Z*(t)] and Z’(t) = g[t 0 , t, Z (t)] , 

a a 


then 


8g[to, t, Z (t)l ^ ( 9g[t, tp, Z'(t)I 1 - 1 

az^Ct) | ~ az’(t) J 


and 


-. ^ L * ^ | p ar tial derivative transition matrix ^ 

evaluated with Z 0 = Z’(t). Thus, the final form of the variation-of-paramel 
differential equation is 

i . 

Z 0 =Z’(t) 

where Z is calculated using Z (t) given by the closed-form solution applie 
P a ~ 

to Z'(t). That is, 


Z (t) = g[t, t 0> Z*(t)] 
a 
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so that the right-hand side of the variation-of -parameters differential equation 
can be seen to depend only on [ t, Z* (t) ] . Thus, numerical integration of this 
differential equation will yield Z’(t) at any t > t^. Then the closed-form 

solution can be used to determine Z (t) = git, t$, Z’(t)] at any t > t^. 

• a 

Since Z^ is usually very small, the differential equation for Z' (t) is usually 


integrated much more rapidly than the differential equation for Z (t) . An 

a 

additional benefit of variation-of-parameters integration as explained here is 


that the transition partial derivative matrix 


9Z (t) 
a 

dZ n 


can be approximated 


9Z a (t f ) 

very easily. To see how this is done, suppose a value of — — - — is desired. 

dZ 0 

Then the interval (to, t^) can be divided into as many subintervals (t 0 , tj) , 


(ti* 


(t , t.) as desired. If this division is fine enough, 
n t 


8 2 » (t ‘ ) „ ( 8 Z M \ 

dZ 0 [dzijtjj 


dZ a (t t } azfa) 
3Z (0 az^iQ 

a 


9Z a (t f ) 9Z(t f ) 

9Z (t ) ” 9Z _, (t“J * 

an f 

Then 


9Z (t ) 
a f _ 

' 8Z(t f ) ' 


9Z (tj>) 


dZ(t<) 

9Z 0 

9Z’(t f ) 


dZ 1 JtJ_ 


9Z f ({) 


Note that the quantities 


r 9z(ti) i 


" 9Z (t>) ' 


'3Z(t f )' 

IzW. 

f 

|_8Z f (t2> _ 

I *'* 9 

_9Z*(t f )_ 


are evaluated 
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using the closed-form partial derivative transition matrix ^ ^ with Z 0 

equal Z f (ti), Z ’(t^), ••• , Z f (t f ), respectively. 


If a numerical integration routine with step-size control is used to 


integrate 


dt 


[ Z* (t) ] , then the sub interval boundaries t if 


* 2 » 


• • • , t will 
n 


be selected automatically to assure that across each interval the approximation 
in the partial derivatives is accurate. This is because the requirement of a 
specified accuracy on the integration of the variation-of-parameters differentia] 
equation will not allow Z to become very large in any particular interval of 


time. Then the beg inni ng of a new interval can just be considered as a new set 
of initial conditions. The following diagram might help to visualize this 
approach: 




w w 


The- diagram exaggerates the differences Z'(t t ) and Z 0 , Z’^) and 

Z (t,), ••• , Z'(tJ and Z (t ); but if the subinterval boundaries 
a 1 f an 

(tj, fc>, • • • , t^) are selected by a numerical integration routine with a 
step-size control (based on accuracy requirements) , the differences will be 


small and the approximations 


az a (t f )' 


azttj) ' 

az (t ) 

a n J 


azuy 


"az (t/ 

a 1 


’az(ti) 1 

’az^tg)' 

. 9Zo 


_az»(ti)J 

’ az (t t ) 

a . 


will be accurate. 




9Z(t,) ' 

az»(t2>J’ 


This completes the explanation of the derivation of a variation-of- 
parameters numerical integration technique to obtain Z (t) and the partial 

az a (t) 3 

derivative transition matrix — — — . In the next section, this information 

dZo 

will be used to solve what is normally called the optimization boundary value 
problem. An algorithm that can solve the optimization boundary value problem 
rapidly enough is, in effect, a guidance scheme. 
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GUIDANCE SCHEME DEVELOPMENT 

The problem of guiding a space vehicle consists of determining an 
optimal trajectory from a measured set of initial state conditions to a specified 
set of boundary conditions. To formulate this problem more precisely, define 
the state vector denoted by X(t) as follows: 


X(t) = 


x(t) 

v(t) 

q(t) 


where to ^ t ^ t 

Then X(t) can be given by the closed-form solutions or obtained from 
numerical integration of the actual nonlinear differential equations. Also 
define the multiplier vector P(tJ as follows: 


P(t) = 


X(t)' 

u(t) 

■y(t) 


Closed -form solutions for X(tJ and P(t) can be used or the non- 
linear differential equations for X and P can be integrated numerically by 
the variation -of -parameters technique to yield X(t) and P(t) . In 


either case the matrices 

r 3x(t)' 

and 

9P(t)T 

L 9 PoJ 

J 

r 9z (t) 

and 

' 8Z a (t) l 

L az ° - 

0Z O J 

roxw 

and 

'9P(t)] 

3Po 

9Po 

the guidance boundary vs 


MW 1 and r 1 ZMT 


are available because 


iPTJ L^TJ 

are a part of the partial derivative transition matrices 

developed in the previous sections. The matrices 
are used to set up a Newton's method iteration to solve 


of physical boundary conditions to be satisfied by a space vehicle as follows: 


F[X(t f ), t f J = 0 . 

Then the optimization boundary value problem consists of determining an 
optimal trajectory X(t) from a measured set of initial conditions X(to) such 
that X(t^) satisfies the physical boundary conditions given above. In order 
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to determine an optimal trajectory X(t) satisfying the above conditions, the 
vector P(t) must be introduced so that the control vector p shown in a 
previous section, Development of Closed-Form Solutions, can be defined in 
terms of X(t) and P(t). Also, from the necessary conditions of optimization 
theory, the following transversality conditions associated with the physical 
boundary condition F[X (tj), t^] must be satisfied: 

, pfc ,T «‘V T / 3FtX < t £ ) ’ ‘f 1 I 
1- P 'V ‘ " p { 8X^3 J ’ 

T 

where J(t £ ) is the scalar function to be minimized and p is a new vector of 

constant multipliers with the same dimension as the vector function 
F[X(t f ), t f ] 

2. J(t ): + p T F[X(t f ), t^i = 0 . 

Now the physical boundary conditions and the above transversality conditions 
can be combined to form the total set of boundary conditions denoted by 

G[X(t f ), P(t f ), p, t f ] = 0 . 

Note that the dimension of the vector function G is equal to the dimension of 
the vector X(t) plus the dimension of the vector F[X(t £ ), t £ ] plus one. Then 

the optimization boundary value problem consists of determining P (t 0 ) , p, and 
t £ such that the boundary Conditions 

G[x(t f ) , p(t f ), p, y = 0 

are satisfied by an optimal trajectory X(t) originating from a set of measured 
initial conditions X(to). The following is a reiteration of the previous state- 
ment of the optimization boundary value problem. 

1. The term X(tg) is measured by the space vehicle, 

2. Values of P(t 0 ), p, and t £ must be determined so that when 

X(tQ) and P(to) are used as initial conditions for the closed-form solutions 
or the variation-of-parameters integration, the resulting solution X(t) and 
P (t) when evaluated at t f will satisfy the total set of boundary conditions 

G[X(t f ) , P(t £ ), p, t f ] = 0 . 
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To solve the optimization boundary problem as stated by the preceding 
sentences, note that X(t { ) and P(t f ) depend on the initial conditions X(t 0 ) 

and P(to). Thus, a Taylor Series expansion of the vector function G[X(t f ) , 

P(t f ), p, t ] can be made about a set of guessed values (denoted by P(to)*, 

P*, and t f *) for the parameters P(to), p, and t { . That is, assume that the 

initial conditions X (to) , P*(t,>) yield X*(t £ *), P*(t £ *). 

Then 


G[X(t f ), P(t f ), p, t f ] = G[X*(t*), P*(t f *), p*. t f * 

9G[X(t ), P(t f ), p, tj 

+ ' — W — L i Ap(t »> 


9G[X(t f ), P(t £ ), p, t f ] 
r ' 9P _ 


Ap 


(d(G[X(t ), P(t ), P, t 1 ) 

+ | 5 


where the subscript * of the braces means that the entire term in the braces 


is evaluated with * values. The terms 


d(G[X(t f ), P(t { ), p, t f ]) 
dt 


9G[X(t f ) , P(t f ), p, t f ] 
— 


and 


can be evaluated explicitly, but the term 


j 9G[X(t f ), P(t f ), p, t f l | 

| j can only be evaluated using the matrices 


ax(t f ) 
9PT ij 


and 


J* 


3P(t f ) 
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That is, 


J.0G[X(t £ ), P(t f ), p, t f ] | 

J 9G[X(t f ), P(t f ), p, t f | 

*ex(t f )' 

{ ‘ ap^5 ' j „ 

( 8X^1 I* 

epl?) 

i j 


j 9G[X(t f ), P(t f ), p, t<] | 

'9P(t f )‘ 

1 ^ 1 , 

8P K) 

- 


Now a modified Newton’s iteration formula can be developed from the above 
Taylor Series expansion truncated after the first-order terms. To do this, 
note that the desired value for G[X(t £ ), P (t £ ) , p, t £ ] is zero so that 


G[X(t f ) , P(t f ) , p, t f ] = a G[X*(t £ *), P*(t f *), p*, t f *] , 

where 0 ^ a < 1, can be substituted into the truncated Taylor Series expan- 
sion for G[X(t £ ), P(t £ ), p, t £ ] to give corrections AP(t 0 ), A p, and At £ 

which, when added to P*(to), p*, and t*, will produce a value of G[X(t f ), 

P(t £ ), P, t £ ] nearer zero. That is, 


a G[X*(t f *), P*(t f *), p*, t f *] 


G[X*(t f *), P*(t f *), p*, t f *] 


f 3G[X(t.) , P(t ), p, tj 

+ | > AP W 


9G[X(t f ), P(t ), p, t.J 

+ 1 ‘—ST Ll Ap 


or 


AP(t 0 ) 
A p 

At f . 


I d(G[X(t f ), P(t f ), p, t f ] ) 


= (a - l)[E] -1 |G[X*(t*), P*(t*), P*. t*] | , 


} 
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where 


r|8G[X(t f ), P(t { ), p, t f ]| MOtX(t), P(t f ), P, t 1 \ | d(G(X(t { ), P(t { ), p, t { !) 1 

El “ [( 8PpJl (« ’ I T }*’( dt ),] 


AP(t 0 )' 


The above expression for the correction vector I A p 


is used repetitively 


until the boundary conditions G[X(t £ ), P(t ), p, t £ ] are zero to a desired 
tolerance. 

The preceding discussion completes the general explanation of the 
solution of the optimization boundary value problem. 

For illustration purposes, this development will now be applied to a 
specific set of physical boundary values. That is, let 


F[X(t f ), t £ ] = v f T v f - v^ 


[ ^ V f - R d V d COS *d J 

Then F[X(t £ ), t £ ] is a three-dimensional vector and the quantities R^, v^, 

and i> , are constant desired values for the position vector, velocity vector, 
d 

and the angle between the two, respectively. Assume that J = -q^ is to be 
minimized. Then the transversality conditions become 


x' T (t f ) = 


0 , 0 , 0 
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u T (t f ) 


= -p 


0 , 0 , 0 


_T 

2V f 


-T 

X f 


y(t f ) 


X 1 (t f ) x(t f ) + ti 1 (t f ) x(t f ) + q = H(t f ) = 0 


_T 


When the transversality conditions are combined with the physical boundary 
condition, the following form for G[X(t^) , P(t^), p, t^] results: 


G[X(t f ), P(t f ), p, t { ] = 


_T _ . 

x. x, - Rj 
it d 


-T - 2 

v v - 
f f d 


_T _ 

x, V. - R v , cos & - 
if d d d 


2 x f> 


■ v f? p+ \ 


V 


0 

n 

, 2v ( , x f \ 

u 

) 


p + u 


-T - _T ^ 

\ X f + U f V f + ^ 


From the above expression for the vector GlX (t^), P(t^), p, t f ] , 
the expression for the components of the [ E] matrix can be obtained. That is 


f8G(X{t,), P(t f ), p, t f ] J 

8G[X(t ( ), P(t f ), P, t f l 


(8G[X(t { ). P(t f ), p, t f j| 

(dG[X(t { ), P(t { ), p. t f l 

1 **• J. 1 

8^ 

* ' 

1 ■*" !•* 

1 
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where 


„ -T 
2x f 


{ 0 , 0 , 0 ] 


9G[X(t f ) , P(t f ), p, t f ] 




-T 

V f 


P 1^3X3 


P 3^3><3 

\ T 


[ 0 , 0 , 0 ] 


-T 

2v f 


_T 

X f 


P 3 [I1 3X3 


P 2^3X3 


IT 

" U f 
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The preceding work completes the explanation of a Newton’s iteration 
procedure for solving the optimization boundary value problem. Note that this 
Newton r s iteration procedure can be used with either the closed-form solutions 
or the variation-of-parameters solutions. 


CONCLUSIONS 

Closed -form solutions which approximate the optimal motion of a space 
vehicle during powered ascent flight and unpowered reentry flight have been 
developed. Also, the entire system of partial derivative transition matrices 
for these closed-form solutions have been developed. This work allows a 
variation-of-parameters integration technique to be used to obtain more 
accurate representations of the motion of the space vehicle in both powered 
ascent and unpowered reentry. Also, the optimal guidance boundary value 
problem is formulated, and a Newton's method algorithm developed for solving 
the guidance boundary value problem is explained using either the closed-form 
solutions or the more accurate variation-of-parameters solution. 

Peter Leung of Northrop Space Laboratories has checked and devel- 
oped a computer program which evaluates the ascent closed-form solutions. 
Table 1 shows how these results compare with the numerical integration of the 
more accurate equations of motion given in Reference 1. The computer 
programming of the rest of the results developed in this report is still in 
progress. When these computer programs are complete, more numerical 
results will be published. 
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TABLE 1. COMPARISON OF THE CLOSED-FORM SOLUTIONS WITH NUMERICAL INTEGRATION 
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TABLE 1. COMPARISON OF THE CLOSED-FORM SOLUTIONS WITH NUMERICAL INTEGRATION (Concluded) 
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